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SIGNALS & SYSTEMS 


e Signals are changes of any quantity w.r.t. time, space or both. 
e Signals can be represented by functions. 
e Functions are defined as maps or associations between objects (Domain & Range). 
e Thus signals can be represented as maps from domain to range. 
e Commonly used Domain set: 
> Set of Integers 
> Set of real numbers 
e Commonly used Range set: 
> Set of complex numbers 


Dr. K.S. Venkatesh (IIT-Kanpur) Signals & Systems Page 1 of 64 


SIGNALS & SYSTEMS 


e Both Domain and Range can be multidimensional. 

e Signals where the domains are of higher dimension (greater than 1), are called multidimensional signals. 
e Signals where the range are of higher dimension (greater than 1), are called multivariate signals. 

e Range set contains all possible values that the function can assume. 

e Complex numbers comprises Range set. 


Theory of Complex numbers 


e Rectangular representation: 
c=at+jb : j=v-1 
a = Real [c]& b = Imaginary [c] 
e Polar representation: 


c —rz0 ir = Va? +b? 


a-rsinü;b-rcos0& 0 = tan !(2) 
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Operations on Complex numbers: 


e Addition: 
ctc =(ata)+j(bt+b) 
e Subtraction: 
c+(-c)=(a-—a)+j(b—-b) 
e Multiplication: 


c.c =(ab—ab)+j(ab+ab ) : rectangular representation 
c.c —r.r 2019 : polar representation 


e Division: 
e Raising to a power: 


e Extracting n^" root of complex number: 
> MWcwill haven dif ferent n'^roots of c 


0 -(n—1)2 
Nc = Vrzo, ,Nrzə,,... NTZQg , ...... ... NrZzQy : Qn ..-. ue14N 
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e System is an entity that takes a signal and converts into another signal. 


€ Asystem is continuous time if its input-output signal is continuous. 


x(t) 


€ Asystem is discrete time if its input-output signal is discrete. 


x[n] 


yIn] 
pL 


e Continuous time system processes continuous time signals x(t) and produces continuous-time outputs y(t). 
e Discrete time system processes discrete time signals x[n] and produces discrete-time outputs y[n]. 
e Different systems differ in the manner in which they relate the inputs and outputs. 
ə Asystem is also a function as it maps input to output. 
e Amap is Injective (one-one) if it takes different input signals to different output signals. 
If xı # x; then yı # y; 
e A map is Surjective (onto) if every signal in the range set will be found as the output for some input signal. 
ə A system can be either Injective or Surjective, or both or neither. 
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ə Asignal must be defined for every value of the argument. 
e Signals that have finite value at every instant of time are finite valued signals. 
e Signals that have values within a certain limit are called bounded signals while signals for which bounds cannot be 
found are called unbounded signals. 
e All physical signals are inherently bounded because any physical source of signals can produce only limited signal 
power. 
e Symmetries are intimately related to signal representations. 
e Signals can have following symmetry properties: 
> Periodicity: x(t) = x(t — T) 
> Evensymmetry: x(t) = x(-t) 
> Odd symmetry: x(t) = —x(-t) 
e Asignal can be decomposed into even and odd parts : x(t) = x,(t) + x,(t) 
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e uj(t) Family : 
» Unit step ua4(t) = u(t) = |n i - : 


> Unit ramp u_2(t) = r(t) = in : : » 0 


e Exponential Family:e”” (s = o + jw) 
> Real exponential signals xit)” 
> Complex (imaginary) exponential signals x(t) = e/ot 


> General complex exponential signals — x(t) = ete 
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e u, [n] Family: 


» Unitstep u. 4[n] = u[n] = ( :n«0 


> Unitramp = u5[n| = rin) = : n<0 


e Discrete Exponential Family: a” (a = re/®) 
> Real exponential signals 
> Complex (imaginary) exponential signals x[n] = e/9^ 
> General complex exponential signals — x[n] — r^. e/2” 


xinler” :U ep 
>r=18&0>0 
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e Asystem is viewed as transformation or mapping of input to output. 
Transformation on Time: 


e On input signal x(t) (defined on time), application of function on time and then a rearrangement of time axis. 
e Examples of transformation: 

> Translation: f(t) =t+T 

» Flip: f(t) = —f(t) 


> Scaling: f (t) = at, wherea is a constant for a #0 


Transformation on Range: 


e Transformation of a signal in the range. 
e Types of transformation: 
> Translation (adding a constant): f (t) = f(t) +a, where a is a constant 
» Flip: f(t) = —f(t) 
> Scaling: f(t) — af(t),a # 0 
e Range and Time (Domain) can be combined. 
e Interchange of range and domain components of a complete transformation does not affect the end result. 
e Change in internal order of the domain or range transformation, changes the end result. 
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ə A system is a mathematical model of a physical process that relates input (excitation) to output (response). 


Properties of system: 


e Memory or Memory less system 
> System is memoryless if its output y(t) ) depends only on its input x(to) at that same instant to. 
> Attenuators and amplifiers are examples of memoryless systems. 
e Causality 
> System is causal if its output y(t) at to only depends upon input x(t): t € t_0 
> System is anti-causal if its output y(t) at ty depends only upon input x(t) : t > to 
> System is non-causal or acausal if y(t) at to depends only upon x(t) : —oo € t — to € oo 


e All memoryless systems are causal but not vice versa. 
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e Stability 
> Asystem is BIBO (Bounded Input/Bounded Output) stable if any bounded input x(t) results in a bounded 
output y(t) 


Ix(t)) € B, : —oo € t € oo zə İy(E)l SB, : —oo € t € oo 
B, and B, are finite and real. 


e Linearity 
> Asystem is linear if it satisfies Homogeneity and Additivity. 
» Homogeneity: 
x(t) ə y(t) > kx(t) > ky(t) : k is complex 


> Additivity: x(t) > y(t) & x (t) y (t) 
x(t) kx (t) 2 y(t) * y (6) 
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e Property of linearity if possessed by a system, is equivalent to a commutation of the linear operations with the system 
itself. 
e Commutation diagram as a test for Linearity: 


x(t) 
LINEAR 
OPERATION 
x(t) 


LINEAR 


OPERATION 


e Commutation diagram as a test for Time Invariance: 


i TIME SHIFT : SYSTEM yü) 
77 Ba 


x(t) y(t) 


SYSTEM TIME SHIFT 


e Every system for which outputs are invariant to the order in which linear (or time shift) operation and system are 
applied are said to be linear (or Time Invariant). 
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e AnLTI system is both linear and time invariant. 
e Representation of LTI system 


y(t) = | x(t)h(t —1)dv : Continuous — time system 


ó(t) h(t) 


——— LTI SYSTEM — 
x(t) E.ə y(t) = x(t) * h(t) 


y[n] = 2. x[k]h[n — kl : Discrete — time system 


k ——oo 


ó[n] h[n] 


— LTI SYSTEM | —” 
x la m 
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e Impulse representation of an input x[n]: 


oo 


Xn] > x[k]ó [n — kl 


k=—0o 
Convolution: 


e For Discrete-time signals x[n] and h[n], their convolution is given as: 


co 


“El. 2 x[k]h[n — kl 


k=—0o 


e The ouput of a discrete LTI system is the convolution of input x[n] with the impulse response h[n]. 


Properties of Convolution: 


ə Commutative: 
h[n] * xinl = x[n] “ hin) 
e Distributes over addition: 
[xy [n] + x2[n]] * h[n] = xi [n] € A[n] + xz [n] * hin) 
e Associative: 
[xx [n] * x2[n]] * h[n] = x3 [n] * [xz [n] * h[n]] 
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e Impulse representation of a continuous LTI system for input x[n]: 


x(t) = b x(1)ó(t — v)dr 


ó(t) h(t) 
area Bu rea SYSTEM 
x(t) 777 y(t) = x(t) * h(E) 


e For Continuous-time signals x(t) and h(t), their convolution is given as: 
y(t) = x(t) s h(t) = | x(t)h(t — t)dt 


e The ouput of a discrete LTI system is convolution of input x[n] with the impulse response h[n]. 
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Properties of Convolution: 


e Commutative: 
x(t) s h(t) = h(t) * x(t) 
e Distributes over addition: 


[x1 (€) + x2(6)] * h(t) = x1 (©) * h(t) x2 (t) * h(t) 
e Associative: 


[xi (€) * x2 (t)] * ACL) = x4 (6) * [x2 (t) * A(t)] 
Evaluation of system properties from the impulse response: 


e Memory less: 

h(t) — kó(t) 
e Causal: 

h(t) 20 :t«0 
e Stable: 


[. |h(t)|dt € oo 
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INTRODUCTION: 


e  Continuous-time, linear, time invariant, dynamic systems are described by linear ordinary differential equation with 
constant coefficients. 
e Mathematical models of such systems: 


x(t): forcing function or system input 
y(t): solution or system output 


e The solution of a linear differential equation y(t) has two additive components: 
> Homogeneous Solution y, (t) 
> Non Homogeneous Solution yı (t) 


y(t) > ys (t) t yas (t) 
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e The homogeneous solution is contributed by the auxiliary conditions of the system. 
e Ahomogeneous solution satisfies the homogeneous differential equation: 


d'y) _ 
Aq 


k=0 


e Solving homogeneous equation: The space of homogeneous solutions has a dimension equal to the order of the 
differential equation. 

e The Non-homogeneous solution is contributed by the forcing function of the system. 

e ANon-homogeneous solution satisfies the Non-homogeneous differential equation: 


N 


ak 
>, Er =O 


k=0 


e Solving Non-homogeneous Equation: This follows from the recognition of classes of functions that remain invariant to 
differentiation and additive combination. Functions such as constants, polynomials and exponentials. Other arbitrary 
forcing functions need to be expressed in terms of these invariant functions. 
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e İn addition to the external forcing function, the system is also driven by its internal forces. 
e This internal force comes from the system's initial / auxiliary conditions. 
e Multiple solution trajectories result from the mathematical solution of differential equations of the system. 
e İn order to identify from the multiple solutions, the unique solution that corresponds to the actual physical system: 
> The description provided by the differential equation is used. 
> The additional information outside the scope of differential equation is required to be supplied. 
» The additional information can be supplied in a variety of forms at instants of time of our choosing. 
> If the additional information is adequate, then a unique solution will be deduced, that matches with the 
behavior of the system. 
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Linearity: 


e Asystem is said to be linear if it meets the following two criteria: 
» Homogeneity: kx(t) > ky(t) 
> Additive: x1(t) ə y4(t) & x;(t) > y2(t) 
x1(t) + x2(t) > y1(0 + yz(6) 


The linear behavior of system described by a differential equation: 


e The response of a system is linear against input if all the auxiliary conditions are zero at the reference instant. 


e The response of a system is linear w.r.t. any one auxiliary condition when input and all the other auxiliary conditions 
are set to zero. 


e Since response consists of two components: 
> One component is linear w.r.t. Input : y; (t) 
> Other component is similarly linear w.r.t. to auxiliary condition: y4(t) 


y(t) > ya) + yi(t) 
: y(t)= Output/response 


: y4(t) =Response when input is zero or Zero input response 
eyi (Öz Response when auxiliary condition is zero or Zero state response 
e The complete response is neither linear w.r.t. to auxiliary condition or Input if the other is also present.. 
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Time Invariance: 


x(t) > y(t) 


or all iç , oo «ty € oo 
x(t — to) > y(t — to) lo all x(t) 


The Time Invariance behavior of system described by a differential equation: 


ə A system described by a differential equation is time invariant when the time of specification of the auxiliary 
conditions is shifted by the same amount by which we shift the input signal. 

e Alternatively, it is time invariant if, instead of shifting the auxiliary condition(s), we modify their values to those that 
they would assume at the new shifted position when we shift the input. 


Causality: 


e The differential equation itself does not compel the system it describes to be either causal or anti causal. 
e The Differential Equation is neutral with respect to the direction of time and can be solved in either direction as we 


choose. 
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e İt isthe counterpart of the differential equation for describing discrete-time systems. 
e The general form of the linear constant coefficient difference equation: 


n 


2, aeyln = K] = xin) 


k=0 
To solve a Linear constant coefficient ordinary difference equation: 
The solution of a linear difference equation y[n] has two additive components: 


> Homogeneous Solutiony, [n] 
> Non Homogeneous Solution y, [n] 


yinl > yalnl £ yaa [n] 
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e Ahomogeneous difference equation: Y ga,y[n — kl 2 0 
e Solving the homogeneous equation: 
> Substituting y[n] for ar”. 
> Factor out the exponential ar”, so one is left with the algebraic characteristic equation. 
> Solve the algebraic characteristic equation and find the roots. 
> For N distinct roots, the solution takes the form 


e Anon-homogeneous difference equation: Y ren ay y[n — kl = x[n] 
e Solving Non-homogeneous Equation: 
> Forcing functions for which we attempt to solve for non-homogeneous equation should be 
constant/exponential /polynomial. 
> The solution of non-homogeneous equation for the above kinds of forcing function will proceed by assuming 
that the solution y [n] is a sequence of the same type. 
> The basic solution of the non-homogeneous equation is denoted by y[n] 
e Thus the complete solution takes the form: 
y[n] = vo[n] + X, cr”: cy :arbitrary constant 
e The relationships pertaining to the specification and insertion of the auxiliary conditions to obtain a unique physical 
solution follow closely those discussed already for the case of the differential equations. 
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Analogous to the continuous time case, the properties of a system can be developed based on approach of differential 
equation. 


Linearity: 
e 


The Response of a system is linear against input if the auxiliary condition is zero at the reference instant. 


e The Response of a system is linear w.r.t. auxiliary condition when input is set to zero. 


y[n] > yalnl + yi[n] 


: y(t) = Output/response 
: y4(t) =Response when input is zero or zero input response 


:y;(t) =Response when auxiliary condition is zero or zero state response 


e Combination is neither linear w.r.t. to auxiliary condition or Input 
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Incremental Linearity 


ə İn the presence of nonzero auxiliary conditions, the output is not linear w.r.t the input as mentioned earlier for 
differential equations. But the system exhibits ‘incremental’ linearity: increments of the output are linear with respect 
to increments of the input. This is because of the cancellation of the output component due to the auxiliary 


conditions. 
x(t) 2 yi(6) : yi(t) € y4(0 +yn® 
x2(t) > y2(t) : y2(t) = WO + yr; (t) 
« [x2(t) —x1(€)] ^ « [yp (©) — yn (Öl 


Time Invariance 
xin — ng] > yin — no] 


e A system described by a difference equation is time invariant when the time of specifications of the auxiliary 
conditions are shifted by the same amount by which we shift the input signal. 
e Alternatively it is time invariant if shifting of the auxiliary condition values that we expect to find at those same 


time(after shifting is done) instants, is done. 


Causality 


e The difference equation is neutral with respect to the direction of time, can be solved both causally or anticausally. 


e Real physical systems are always causal. 
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Moving Average Filter 


e It essentially creates a gradually moving average, at every instant of time the window over which the input sequence 
is being averaged is being shifted by one position, along with the point of interest. 
e It is called a Finite Impulse Response (FIR) filter as (the support of) its impulse response is finite. 
e The order of filter depends on the number of delay blocks. 
e The general form of an M“" order moving filter : 2 b, x[n — kl) 
e The finite impulse response of Moving Average filter : h[n] = X4 o b, öln — kl 
Auto Regressive Filters: 
ə For Auto Regression filters of order N, the system output at time n is a linear combination of N past outputs. 


e It is called Infinite Impulse Response (IIR) filter as its impulse response is invariably (of) infinite (support). 
e The general form of an Auto Regressive filter : Xo a, y[n — kl 


Auto Regression Moving Average (ARMA) filters: 


e İt is combination of Auto Regression (IIR) & Moving Average (FIR) filters. 
e The general form of IIR & FIR filter combination: Y oa, yin — kl = Xİ” b, xin — k] 


Filter Implementation: 
> Direct Form 1 
> Direct form II 
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e Integrators are an alternative to differentiators for continuous time implementation of the solver of differential 


equation. 
Integrators average out minor fluctuations or noise unlike differentiator which enhances them with repeated 


differentiation. 


Solving the differential Equation using Integrator: 


N N 
d*y(t) | " d* x(t) 
LET * dik 

k=0 KO 


e Equation after Integration over N time : 3 a, yln-*l(t) = Ned b, xln-*l(t) 
e Final equation : V(t) = sb, xin-kl(t) — 2. a, yl'-*l(t) 


Filter Implementation: 


> Direct Form 1 


> Direct form Il 
e Fora causal implementation of a difference/differential equation solver, the auxiliary condition necessarily has to be 


the initial condition. For an anti-causal system being solved for t o —oo to t = 0, the auxiliary condition has 


necessarily to be a final condition. 
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e The motivation for signal representation is the convenience of dealing with a small number of signals taken out of a 
fixed collection. 

e This provides a much more compact description of the systems behavior. 

e Instead of dealing with a processor in terms of what it does to each possible signal in the signal space, it is much more 
convenient to decompose the signal space into small classes. 

e Each of these small classes is preserved under action by the processor in question: the output is also in the same 
class as the input. 

e Arbitrary signals outside of these preserved classes can be represented as a combination of members of the 
preserved classes. 

e The combining process commutes with the processor: eg, additive combining commutes with linear systems. 

e Thus, effectively, applying the processor directly upon the input signal is equivalent to a decomposition under the 
representation, applying the processor upon each component, and combining the results to get the output. 

e The choice of a signal representation scheme is strongly dependent upon the kind of systems we are concerned with. 
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e Acontinuous time signal x(t) is periodic if there exists a positive and non zero T such that 
x(t +T) 2 x(t): for all t 
e Given x(t) 2»x4(t) +x2(t) : xı(t) € x(t —T) & x2(t) = x2(t — T2) 
e Signal x(t) will be periodic if 
T. 1 
» - —nor z such that 71 is an integer and periodicity T = maximum ( T3, T; ) 
2 
T 


+= “then periodicity T = LCM ( T, T; ) 
T2 n 


> 


. : “an ; 
e Signal x(t) remains non periodic if - z rational number 
2 


e Given a period T, another period T' is said to be harmonically related to it when T/T’ is a positive integer. 
e Ifa T-periodic signal is to be represented as a sum then, 
> Allrepresenting components must be periodic. 
» Additionally, all the periods of representing signals must be harmonically related to the period of x(t). 
ə Apreserved representation class for T-periodic signals in the context of processing by LTI system, will be of the form 
(e tH, k = 1,2,3...} 
> The signal is a complex exponential and its form is preserved by any LTI system. 
> It is periodic and harmonically related to that of x(t). 
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e İt is a linear representation for periodic signals that are real or complex valued. 
e Fourier series representation of a periodic signal x(t) (Synthesis): 


oo 


.2T 
xit —x(6—T)— 2. x,e T^ 
k-——oo 
e Fourier co-efficient x, (Analysis): 
‘ y 
əm 
= JEO oir di 
0 


Convergence (in energy) of Fourier Series: 
7 T 

If fy lx (t)|?dt € co : x(t) has finite power 

=> all coefficients xı will exist. 


— difference between x(t) and its representation will have zero energy. 
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Convergence (pointwise) of Fourier series: 


e Dirichlet conditions: If x(t) is continuous at a point, the FS representation will converge 
e Periodic signal x(t) has a Fourier representation if it satisfies following Dirichlet conditions: 
> x(t) is absolutely integrable over a period 


T 
| |x (t)|dt « oo 
0 


> x(t) has at most a finite number of maxima and minima over one period. 
> x(t) has at most a finite number of discontinuities over a period and each discontinuity is finite. 
e Atadiscontinuity, the FS representation converges to the midpoint of the discontinuity. 


Properties of Fourier series: 


e Fourier series is an orthogonal representation. 
e Fourier series coefficients are generally complex numbers. 
e Ifx, — x', , then x(t) is real. 


e If x, = —x*, , then x(t) is imaginary. 
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e Fourier spectrum is plot of x,, (complex amplitude) for each k*^ harmonic. 
e The Fourier spectrum is often called a line spectrum because at non-harmonic frequencies, the energy is zero. 
e The lowest spectral line is located at @ = « which is the fundamental frequency of the spectrum. 


2n : : . : 523 
ə dm, where T is period of a continuous time periodic signal. 


e Aspectrum can be finely resolved by changing the periodicity of the signal without changing its shape, padding it with 
extra zero space and thus increasing the period. 


e For spectrum given below, Let x(t) has period T & x'(t) has period 4T 


e Spectrum of x'(t) will be 4 times more finely resolved than that of x(t) because its period T” is 4 times T, the period 
of x(t) 
ə Asperiod T > oo, a continuous time periodic signal becomes non periodic and its spectrum becomes continuous. 
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e Fourier transform is a Fourier representation of non-periodic signal x(t). 


e Analysis equation : 


X(w) = (aoa x(t)heJ® dt : forward transform 


e Synthesis equation: 
x(t) = T ff” X (əefet do : inverse transform 


(0 ——00 


Continuous Time Fourier Transform (CTFT) 


e It finds representations for non-periodic continuous time signals, in terms of complex exponentials that are 


themselves periodic but not harmonically related. 
e The number of representing signals e/?* is uncountably infinite. Hence the linear combination that is the synthesis 


equation takes the form of an integral, rather than a sum. 
e The representation is thus in the form of a continuous (not a line) spectrum and is generally complex valued (just like 


Fourier series) and has a real and imaginary components or magnitude and phase. 
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Convergence of Continuous Time Fourier Transform: 


e Any continuous signal x(t) which has finite energy can be represented in frequency domain via Fourier Transform 
which is known to converge in energy to the signal. 


T 
İr (Ol”dt « oo 
0 


e The condition for pointwise convergence of X(w) is the set of Dirichlet conditions: 
» x(t) is absolutely integrable over a period 


T 


İ |x (t)|dt « oo 


0 


> x(t) can have at most a finite number of maxima and minima over any finite interval. 

> x(t) can have at most a finite number of discontinuities over a period and each discontinuity must be finite. 
e Fourier Transform does not exist for finite power signals. It exists for finite energy signals only. 
e Finite Power signal in time can be represented by Fourier transform if we allow impulses in the frequency domain. 


Dr. K.S. Venkatesh (IIT-Kanpur) Signals & Systems Page 33 of 64 


SIGNALS & SYSTEMS 


e The continuous-time Fourier transform may be viewed as a one-one relationship between complex valued functions 
of continuous time and complex valued functions of continuous frequency. 


» Linearity: 
x(t) o X(o) 
ax(t) + by(t) o aX(o) + bY(o) 
» Conjugation: 
x'(t) o X'(—o) , x*(t) : complex conjugate of x(t) 

» Time Reversal: 

x(-t) o X(—o) 
» Frequency Reversal: 

X(—o) e x(-t) 
» Time Scaling: 


x(kt) o at (2) 


Dr. K.S. Venkatesh (IIT-Kanpur) Signals & Systems Page 34 of 64 


SIGNALS & SYSTEMS 


» 


Frequency Scaling: 

X (2) e |k|x(kt) 

k 

Symmetry: 
If x(t) is real : X(w) is conjugate symmetric 
If x(t) is imaginary: X(w)is conjugate anti symmetric 
If x (t) is real conjugate symmetric x(t)is real & even: X(w)is even & real 
If x(t) is imaginary & conjugate antisymmetric: X(w) is conjugate antisymmetric & imaginary 


Time Shift: 
x(t — ty) o e Jot X(o) 
Frequency Shift: 
X(w — Wo) e el “9 x(t) 
Convolution Property: 
H(w).X(w) © x(t) * y(t) 
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» Duality: 
FT[X(t)] = 2nx(—o) 
> Modulation property: 


1 
x(t). ¥(0) & ==X(w) * Y(o) 


> Parseval Relation: 


[. lxCOI dt = x) meo də 


—oo 


> Differentiation: 


dx(t) 
= jwX 
£ 7? (o) 
» To determine running integral of any x(t) 


ras x(t)dt = x(t) s u(t) ; u(t) is unit step 


X(o) 
Y(w) = X(w).U(w) = "xu OO) ; 16(w) =0;w #0 
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e The output of an LTI system with the impulse response h(t) due to input x(t) is given by the convolution. 


y(t) = x(t) * h(t) 
Y(w) = X(w). H(w) 


Frequency Response H(w) = ə 
LTI System Characterized by Differential Equation: 
x . (t) = yə E (0 : MEN 
7 kgqee o 505 
k=0 k=0 
e Taking Fourier transform on both the sides 
N M 
2a. Go)'¥(@) = 9 by Go)  xta) 
k=0 k=0 


¥(w) Theo by Qeyr 
X(w) ” Shaya Gaye WO 


e H(w) is in polynomial rational form so it can be expanded into partial fractions. 
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Cascade form Implementation: 
Factorizing Y (c) and X(w): if the system is real, 


N : "v 
b 1+b4,jw+b o) ə 
H(w) = pal iz Tdbudurbutiay A biz, bək, dik, az, are all real coefficients 


k=1 1-Fd1k/ o bazk (jw)? 
Parallel form Implementation: 


Determination of Partial Fraction Expansion of H (cə) 


N/, 
box + bi jo 
(o) oe): 2 ko) Ko) 1 + ayy jw + Az, (jo)? 


K 
Ho(w) = 2: doe Go! 
k=0 
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e Periodic Discrete Time sequences: x[n] = x[n — N] 


Theory of representation 
e Preserved classes of complex exponentials (discrete) that are harmonically related to x[n] should be considered. 
e Aperiodic complex exponential satisfies: 


,2nk .2nk 
e = el (n-N) 


e The representation of the periodic signal x[n] = x[n — N] can be done with no more than N distinct periodic 
complex discrete exponentials. 


e Synthesis Equation: 
N-1 
jen 
xin] - 2 xem): en] e^ 
k=0 


e Analysis Equation: 


N-1 N-1 
me = p ble" 2 LY xinl oi inl 
N 0 ii 0 
n= n= 
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Matrix form representation of Synthesis and Analysis equation 


Synthesis Equation: x[n] = Xİ” x, e, [n] 
ool) gA[O| — seins Q4 [0] x[0] 
E HK Qo [o eil] ......... aa ; 
-———R 0-10 
Analysis Equation:x, = ua o) [n] 
" 90[0] | eol] ....... PolN — 1] x[0] 
° ej[0] — eil] ......... 95[N — 1] 
iu 93410] es ali] n eis IN — 1] x[N — 1] 
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The DTFT X (Q) of a discrete time signal x[n] is a function of a continuous frequency 2. . 


e Analysis Equation: 


X(Q) = 2 x[n] e 79» 


n 


e Synthesis Equation: 


x[n] = z; | (mea 


—T 


e For YXalxinll € o», X(Q) will be finite for all Q : 0 € Q x 2z 
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e Discrete-time Fourier transform is one to one map of discrete time sequence and periodic continuous frequency 


functions. 
» Linearity: ax; [n] + bx;[n] ə aX4(Q) + bX;(Q) 
» Conjugation: x*[n] e X*(-Q) 
> Time Reversal: x[-n] e X(-Q) 
» Time/ Frequency Shift: xin — ny] e e?9?"ox(Q) 


> Differencing & Summation: 
First difference of x[n] = X(0)(1 — e/®) 


X(Q) 
(=e) 


Running sum of X(Q) : X^, x[n 1 = X(Q). U(Q) = + nX(Q) Xi2-c ö(0 — 21) 


> Differentiation: 


d | 
3g X (0 € x[n](—jn) 
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» Symmetry Properties of the DTFT: 
X(w)is periodic with period 27 
If x(t) is real : X(w) is conjugate symmetric 
If x(t) is imaginary: X(w)is conjugate anti symmetric 
If x (t) is real conjugate symmetric : x(t)is real & even: X(w) is real & even 


If x(t) is imaginary & conjugate anti symmetric: X(w) is imaginary & conjugate anti symmetric 


» Convolution: 


Linear convolution of two periodic signals is not possible. Thus periodic convolution of two periodic signal of same 
period, T: 


N-1 


x[n] 2.5. ]h[n — 1] 


» Modulation Property: 


1 TU 
x[n]. y[n] = zl X(05Y(0 — 05407 


Dr. K.S. Venkatesh (IIT-Kanpur) Signals & Systems Page 43 of 64 


SIGNALS & SYSTEMS 


e The output of an LTI system with the impulse response A [n] due to input x[n] is given by the convolution. 


yinl  x[n] * h[n] 


Y(Q) = X(0). H(Q) 


Frequency Response H (0) = m 
LTI System Characterized by Difference Equation: 
N M 
2, ae yln - kl = 2, bein = 4] :MEN 
k=0 k=0 
e Taking Fourier transform on both the sides 
M N 

pi b, e )*?Xx(Q) = 2. a, e J*9y (Q) 
k=0 k=0 


YO). EB evt. HO) 
X(Q) Yi a, CTO 


e H(0) isin polynomial ration form so it can be expanded to partial fraction. 
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Cascade form Implementation: 
Factorizing Y(Q) and X(Q) 


dz 
b 14+ bie 79 + p,,e 79 
H(Q) =~ 1k 2k 


" 1+a..e7 4 ag, e-2in : biz, bək, Qjy,Q02, are all real coef ficient 
0 k-İ 1k 2k 


Parallel form Implementation: 


Determination of Partial Fraction Expansion of H(w) 


Ni, 


HQ) = Ho(@) +) H( : HO) = 
k=1 


Ho(o) = 2. qo, e 7*9 
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e When the sampling signal (impulse train) d(t) is applied on a band limited signal x(t), the sampled signal x,(t) 


obtained is given as: 
- FI 2 


ttt 
TX t “m T3373 77; After Sampling 
>» 


Impulse train sampled signal 


Input signal 


x(t) = x(t).d(t) = ” x(nT,).8(t —nT,) 


n 


1 
X;(w) = T[x(t).d(t)] = ə [X(w) * D(o)] 
2T 2T 


1 
X9) ==) X(0—ka,) : a — 77 
S k S 
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It is possible to retrieve original signal from its ideal sampled signal using an ideal low pass filter [L(w)] of bandwidth 


e 
B and gain T, such that: 
T, ; ləl x 
L(o) = 
0;|o| z 


Sol S ro | og 


1 
X,(w).L(w) = -2. X(w — kw,).L(w) = X(o) 
Tex 


Sampling theorem states that signal x(t) of bandwidth B can be reconstructed from its samples if sampling is done at 
frequency w, > 2B or T, < " 
If the signal is undersampled, i.e. w, € 2B , the effect of spectral overlapping is called aliasing. When this occurs, the 


the reconstructed signal is distorted with respect to the original. 
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Problems with ideal sampling and retrieval: 
e Impulse trains are the ideal sampling signals, but they cannot be generated practically, so an ideal impulse train 
sampling is not feasible. 
e Similarly, an ideal low pass filter is the perfect reconstruction agent, but is practically impossible to achieve: a real 
(finite order) lowpass filter has transition a band of non-zero width, which picks energy from neighbouring aliases.. 


Solution: 


e Use of rectangular pulse train instead of an impulse train. 
e There are two methods of sampling using rectangular pulse train: 
» Flat-top sampling 
» Faithful-top sampling 
e Use of a hold circuit (finite order interpolator) instead of an ideal low pass filter. 
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Flat top sampling: 


e Rectangular pulses assume a height equal to the value of x(t) at the points of occurrence of the center of the 
rectangular pulses and their top will be flat. 

e Mathematically, a convolution of ideal impulse sampled signal x,(t) with a rectangular pulse will result in flat top 
sampling. 

1 , İt — nT,l ST m. Flat təp 

0 : elsewhere 3 tes pl. rii 

x(nT,) ; İt — nT,l < T 

0  ; elsewhere 


e Rectangular pulse : rr(t) = 


e Faithful top sampled train: x/(£) = b 


xy (t) = x,(t) err(t) 


1 
Xi(o) = 5 X, (o) Rr) 


2 sin wT 
Rr(@) = ——— 


To recover the original signal from the flat top sampled signal: 


e An ideal low pass filter would not recover original signal without distortion. 
e Due to flat top pulses the spectrum of the sampled signal is distorted. 
e The narrower the pulse width (the more rz (t) resembles ó(t)), and hence, the less distortion. 
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e Pulses resulting from faithful-top sampling will follow exactly the values of the input signal during the ON time of each 
pulse in the train. 

e Mathematically, a multiplication of the input signal with a rectangular pulse train is carried out, with the rectangular 
pulse train viewed as a convolution of a single rectangular pulse with an ideal unit impulse train. 

e Faithful top sampled train: X, (t) = x(t). p(t) 


e Rectangular pulse train: p(t) = X, rr(t — nT,) Fas lh] top Sampling 
xi" 
X,() = x().pO =) x(O.rr(t — nT,) 


n 


X.(o) = Fİx(O.p(Ö) = 1X) 4 PO) 


p B 27 omer km) 
2.77” 


1 2sinoT 
X,(0) ==) = X(u — kes) > Ws. = = 
S k o 
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To recover original signal from faithful top sampled signal: 


T, ;|o| € = 
e Filtering it with low pass filter L(w) = K 
0;|o| = - 
X,(w). L(w) = 2X(o) 
e Faithful top sampling followed by filtering it with ideal low pass filter will yield undistorted output 2X(w). 
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e Ideal low pass filters cannot be constructed without investing infinite resources either in the form of hardware or 
computational effort. 

e Problem can be overcome by use of hold circuit which is a non-ideal low pass filter.. 

e Adesirable property of an interpolator system is that the output should coincide with the input at the data points 
given. 

e First order hold circuit acts as a linear interpolator 

e Both zeroth order hold [ho(t)] & first order hold [h4 (t)] circuit uphold the above property of an interpolator because 
h(t) 20; t-nT;;nz0 

e But a second order hold circuit A;(t) onwards this property can be lost. 

e First order circuit produces a series of ramps by connecting tops of these samples by straight lines. 

e Ideal low pass filter acts as a hold circuit of infinite order, but maintains the property of the output of the interpolator 
coinciding with the data at given points by virtue of being self orthogonal for nonzero shifts that are multiples of T,. 
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e An absolutely integrable signal can be represented via its Fourier transform. 

e Many important signals do not have Fourier transform, because the integral that defines X(w) diverges. 

e The remedy can be a Laplace transform of such signals in continuous time. 

e For a continuous time signal x(t) and complex frequency s = (o + jw) the Laplace transform X(s) will be given as: 


oo 


Xi iEoLaL :s—ocjo 


—oo 


e It isa Fourier transform representation of continuous time signal x(t) with a (possibly) decaying exponential 
terme”. 


CTFT 2. 
X(s) — x(t)e € 
e If Laplace transform of a causal x(t) exists for o4, it will also exist for o? where oş > 0 


e The range of values of complex variables s for which the Laplace Transform converges is called the Region of 
Convergence (ROC). 
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ə Laplace Transform, a function of two variables, o, w is mapped on the so-called s plane. 
e Real part of s is represented on the x-axis while Imaginary part of s on y-axis. 
e The Fourier Transform of a signal is the Laplace Transform evaluated upon the jw axis 


Example of s-plane plot with ROC: 


Region of 
Instability 
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The inversion of Laplace transform is determination of x(t) from its Laplace transform X(s). 


» Inversion Formula: 


so misə 
1 
t) - — X SITE 
x(t) öm | (s)e S 
S=o0—jw 
where o is any point in the ROC of X(s) 
> By Partial Fraction expansion: 
If X(s) is rational polynomial function of the form 
X(s) = N(s) _ (s — 21) (S — zş)............. (S — Zp) 
D(s) (s — p3)(s — pə)... vee sees (S — Dy) 


Whens > z,, X(s) > 0;  z,arecalled zeros 
When s > p,, X(s) Ə co, pa are called poles 


The form can be expanded into partial fraction. The simplified sum of X(s) generated can be then be represented as 
inversion representation, where each term inverts to a standard time exponential. 
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Properties of ROC: 


> ROC does not depend on w but only on o. 
Since the ROC is the set of points "s" at which X(s) converges and since X(s) clearly does not converge at s — p, , it 
follows that no poles can lie within an ROC as that would violate the definition of ROC. 


v 


If x(t) has finite support and is absolutely integrable then its ROC is the entire s-plane. 

x(t) is said to be left sided if there exists ty such that x(t) = 0 for t > to. E.g. - Anticausal signals. 

x(t) is said to be right sided if there exists ty such that x(t) = 0 for t € to. E.g. - Causal signals. 

When x(t) is right (left) sided and x(t)e "* is absolutely integrable at s; = op + jwo then the ROC of x(t) extends 
indefinitely to the right (left) of sg and is called a "Rightward (Leftward) ROC". 

When x(t) is two sided then it may be expressed as a sum of left and right sided signal components, and the rationale 


VV V V 


v 


given above applied jointly upon the parts. 
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» Linearity: 
xi(t) > Xi(s) : Ry 
xə(t) > Xə(s) : Ro 
x(t) = x1 (t) t x2(t) 9 X(s) = Xi(s)-X;(s) : RDR, NR 
> Time Shift: 


x(t — ty) o X(s)e 5 :R 


» Frequency Shift: 
X(s—sg)ex(t)e** :R — R koy : 59 — oy 4 jog 


» Time Scaling: 


x(at) e x (*) : R= (E+ jo: o * jo c R] 
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» Convolution: 
xi(t) * xə(t) o Xi(s).Xə(s) : RAR, NR, 


» Differentiation in time: 


S=ot+jw 

1 

x(t) -— İ X(s)e”St ds 

2nj 

S=o0—jw 
dx(t) E inem 
e : 
dt SAS = 


» Running Integral: 
t 


[ sat =u «20 = =X(5) : R =RN (o 1İa,o 5 0) 


—oo 
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e DTFT does not exist if signal x[n] is not absolutely summable. 

The remedy for the problem is Z transform. 

e Ztransform is the counterpart of Laplace Transform for discrete time signals. 

Z transform is to DTFT is what Laplace Transform is to CTFT. 

For a discrete time signal x[n] and complex value z = re/® ; r € [0, co] the Z transform X(z) is given as: 


X@= >, tele ; z — rel) for re [0,0] & Ne [0,27] 


X(z) —— x[n]r ^? 
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e ZTransform, a function of two variables, r, Q is mapped on the so-called z plane. 
e Magnitude part of z is represented on the r-axis while phase part of z on O-axis. 
e The discrete Fourier Transform of a signal is the Z Transform evaluated upon |r| = 1, the so-called unit circle. 


Example of z-plane plot with ROC: 


Im x(e”) 


“AN, Unit Circle 
- D t 


e The set of values of z for which the z-transform converges is its Region of Convergence (ROC) 
e Each value of r represents a circle of radius r. 
e The region of convergence is made of circles. 
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e The inversion of the Z transform is the determination of x[n] from its Z transform X(z). 


» Inversion Formula: 


1 rel" 
x[n] = m) | X(z)z"-! dz 
ren 
» Partial Fraction expansion: 
If X(Z) is of polynomial rational form 
X(s)= N(z) _ (z — z)(z — zə)... .......... (Z — Zn) 
D(z) (z — p))(Z — p) wee see see .. (Z — ps) 


Carry out a partial fraction expansion to obtain smaller terms that can be inverted from standard expressions. 
Assign to each component term an ROC such that the combined ROC of the component terms matches the 


ROC given for the overall expression. 


Dr. K.S. Venkatesh (IIT-Kanpur) Signals & Systems Page 61 of 64 


SIGNALS & SYSTEMS 


Properties of ROC: 


» 
» 
» 


The ROC of the Z Transform consists of complete circles about the origin. 

The ROC cannot contain any poles. 

If x [n] has finite support İn,, n; ]. then x[n] is always summable and so will be x[n]r ^". If the support of 

x[n] includes points where n is negative, then z = oo will not lie in the ROC. If the support of x[n] includes points 
where n is positive, then the origin will be a point of nonconvergence. 

If x[n] is right (left) sided and the ROC of X(z) contains the zy = ry e/?^, then it will also include all points for z 
satisfying |z|»|zo| (|z|»|zol). Further, if the support of x[n] includes points where n is negative, then z = oo will 
not lie in the ROC. If the support of x[n] includes points where n is positive, then the origin will be a point of 
nonconvergence. 

Two sided signals may be expressed as a sum of left sided and right sided signals. The ROC will be the intersection 
of the component ROCs if they both exist and intersect. 
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» 


Linearity: 

ax,[n] + bx,[n] ə aX,(z) + bX,(z) : R2hR,nR, 
Time Shift: 

x[n — ng] e z-"°X(z) : R) = R except possibly for z = 0 
Frequency Shift: 


X[z — zo] € z"x[n] : R) =Rexcept forz — 0 


Time Reversal: 
x[n] ə X(z) : R 


x[^n] + x (2) GR =; 


Convolution: 
xi[n] * x2[n] > X1(2).X2(2 : RAR NR, 


Differentiation in the Z-domain: 


ax 
ə = > x[njz"1(-n) =: R 
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